issues: a code structure that allows low 
encoding complexity, and a flexible rate- 
compatible code that allows matching to 
various channel conditions. Most of the 
previous high-performance LDPC codes 
for the relay channel are tightly opti- 
mized for a given channel quality, and 
are not easily adapted without extensive 
re-optimization for various channel con- 
ditions. This code for the relay channel 
combines structured design and easy en- 
coding with rate compatibility to allow 
adaptation to the three links involved in 
the relay channel, and furthermore of- 
fers very good performance. The pro- 
posed code is constructed by synthesiz- 
ing a bilayer structure with a protograph. 
In addition to the contribution to relay 
encoding, an improved family of proto- 
graph codes was produced for the point- 
to-point AWGN (additive white Gaussian 
noise) channel whose high-rate mem- 


bers enjoy thresholds that are within 0.07 
dB of capacity. 

These LDPC relay codes address three 
important issues in an integrative man- 
ner: low encoding complexity, modular 
structure allowing for easy design, and 
rate compatibility so that the code can 
be easily matched to a variety of channel 
conditions without extensive re-opti- 
mization. The main problem of half-du- 
plex relay coding can be reduced to the 
simultaneous design of two codes at two 
rates and two SNRs (signal-to-noise ra- 
tios), such that one is a subset of the 
other. This problem can be addressed by 
forceful optimization, but a clever 
method of addressing this problem is via 
the bilayer lengthened (BL) LDPC 
structure. This method uses a bilayer 
Tanner graph to make the two codes 
while using a concept of “parity forward- 
ing” with subsequent successive decod- 


ing that removes the need to directly ad- 
dress the issue of uneven SNRs among 
the symbols of a given codeword. This 
method is attractive in that it addresses 
some of the main issues in the design of 
relay codes, but it does not by itself give 
rise to highly structured codes with sim- 
ple encoding, nor does it give rate-com- 
patible codes. The main contribution of 
this work is to construct a class of codes 
that simultaneously possess a bilayer par- 
ity-forwarding mechanism, while also 
benefiting from the properties of proto- 
graph codes having an easy encoding, a 
modular design, and being a rate-com- 
patible code. 

This work was done by Dariush Divsalar of 
Caltech, and Thuy Van Nguyen and Aria Nos- 
ratinia of the University of Texas at Dallas for 
NASA’s Jet Propulsion Laboratory. For more in- 
formation, contact iaoffice@jpl.nasa.gov. 
NPO-47539 
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For numerical simulations of such flows, fine-grid LES is not as accurate as coarse-grid LES. 
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Multiphase turbulent flows are en- 
countered in many practical applica- 
tions including turbine engines or natu- 
ral phenomena involving particle 
dispersion. Numerical computations of 
multiphase turbulent flows are impor- 
tant because they provide a cheaper al- 
ternative to performing experiments 
during an engine design process or be- 
cause they can provide predictions of 
pollutant dispersion, etc. Two-phase 
flows contain millions and sometimes 
billions of particles. For flows with volu- 
metrically dilute particle loading, the 
most accurate method of numerically 
simulating the flow is based on direct 
numerical simulation (DNS) of the gov- 
erning equations in which all scales of 
the flow including the small scales that 
are responsible for the overwhelming 
amount of dissipation are resolved. 
DNS, however, requires high computa- 
tional cost and cannot be used in engi- 
neering design applications where itera- 
tions among several design conditions 
are necessary. Because of high computa- 
tional cost, numerical simulations of 
such flows cannot track all these drops. 

The objective of this work is to quan- 
tify the influence of the number of com- 
putational drops and grid spacing on 


the accuracy of predicted flow statistics, 
and to possibly identify the minimum 
number, or, if not possible, the optimal 
number of computational drops that 
provide minimal error in flow predic- 
tion. For this purpose, several Large 
Eddy Simulation (LES) of a mixing 
layer with evaporating drops have been 
performed by using coarse, medium, 
and fine grid spacings and computa- 
tional drops, rather than physical drops. 
To define computational drops, an inte- 
ger Nr is introduced that represents the 
ratio of the number of existing physical 
drops to the desired number of compu- 
tational drops; for example, if ,Yr=8, 
this means that a computational drop 
represents 8 physical drops in the flow 
field. The desired number of computa- 
tional drops is determined by the avail- 
able computational resources; the 
larger Nr is, the less computationally in- 
tensive is the simulation. A set of first 
order and second order flow statistics, 
and of drop statistics are extracted from 
LES predictions and are compared to 
results obtained by filtering a DNS data- 
base. First order statistics such as Favre 
averaged stream-wise velocity, Favre av- 
eraged vapor mass fraction, and the 
drop stream-wise velocity, are predicted 


accurately independent of the number 
of computational drops and grid spac- 
ing. Second order flow statistics depend 
both on the number of computational 
drops and on grid spacing. The scalar 
variance and turbulent vapor flux are 
predicted accurately by the fine mesh 
LES only when Nr is less than 32, and by 
the coarse mesh LES reasonably accu- 
rately for all .Yr values. This is attributed 
to the fact that when the grid spacing is 
coarsened, the number of drops in a 
computational cell must not be signifi- 
cantly lower than that in the DNS. 

Results indicate that for large .Yr val- 
ues, the fine-grid LES is not as accurate 
as coarse-grid LES, besides being com- 
putationally more intensive. This is at- 
tributed to the fact that a fine-grid LES 
used in conjunction with a reduction in 
the number of followed drops from the 
physical to a computational drop field 
implies that there is necessarily a smaller 
number of drops in a computational cell 
than in DNS; this aspect naturally influ- 
ences the flow development and biases it 
from the filtered DNS. The key to suc- 
cess seems to be having approximately 
the same and not a significantly smaller 
number of drops in the computational 
cell volume in LES as compared to those 
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in DNS. Thus, the notion that smooth 
statistics are all that is required in the 
choice of the number of computational 
drops does not hold. This is because un- 
like particles in a Monte-Carlo calcula- 
tion where the number of these tracked 
particles would be increased to obtain 


convergence of results, drops are not 
stochastic particles as they obey a deter- 
ministic trajectory that is only modu- 
lated by turbulence. This finding also 
highlights the importance of obtaining 
drop-number experimental data in com- 
bustion chambers or natural flows. 


This toork was clone by Josette Bellan and 
Senthilkumaran Radhakrishnan of Caltech for 
NASA’s Jet Propulsion Laboratory. For more in- 
formation, contact iaoffice@jpl.nasa.gov. 
NPO-48202 
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